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Abstract 



Direct photon production in longitudinally polarised proton-proton 
collisions offers the most direct and unproblematic possibility to de- 
termine the polarised gluon distribution of a proton. This information 
could play a major role for improving our understanding of the nucleon 
structure and QCD in general. It is hoped that such experiments will 
be done at RHIC. We present results of detailed Monte Carlo simu- 
lations using a code called SPHINX. We find that for RHIC energies 
and large gluon polarisation the Compton graph dominates allowing 
for a direct test of /S.g. Triggering on away-side jets with the envis- 
aged jet-criteria should allow to obtain more detailed information on 
^g{x). The photon asymmetry resulting from the asymmetry of pro- 
duced provides an additional signal, which is complementary to 
the other two. For small gluon polarisation, i.e. A^f < 0.5 or very soft 
polarised gluon-distributions the envisaged experiments will require 
a highly sophisticated simulation and large statistics to extract more 
than upper bounds for |Ag'(x)|. 



PACS-No.: 13.88.+e, 13.85.Hd, 12.38.Bx 
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1 Introduction 



The experiments on longitudinally polarised lepton-nucleon scatter- 
ing performed over the last years at CERN j^, |2[ and SLAC have 
determined the polarised structure function gi{x) with rather good 
accuracy, while gi{x) still needs improvement. The theoretical inter- 
pretation of these data is the subject of intense debates which we do 
not want to review here. These discussions have, however, established 
a number of facts which imply the importance of analysing polarised 
proton-proton collisions: 

1. ) Polarized reactions provide very sensitive tests of QCD. To re- 
alize this potential fully requires, however, a detailed experimental 
knowledge and theoretical understanding of Q^-dependences, just as 
for unpolarised reactions Q. To disentangle the various effects pre- 
cise and complementary measurements will be needed. In addition 
comparing specific results from deep-inelastic scattering and proton- 
proton collisions allows for interesting tests. 

2. ) The polarised gluon distribution plays a very special role as it con- 
tributes via the anomaly Q in the same way as the polarised quark 
distribution. For deep inelastic scattering there is no way to distin- 
guish between a 'genuine' quark distribution and an 'anomalous' gluon 
distribution. Such a distinction is by principle only possible by com- 
bining inclusive and exclusive data, see e. g. for polarised deep- 
inelastic scattering, though probably very difficult in praxis. On the 
other hand to establish the existence of an anomalous gluon contri- 
bution would test fundamental topological properties of QCD, which 
have never been directly accessible before. 

Obviously in this situation a precise knowledge of /S.g{x^Q'^) would 
help tremendously, and this could be provided by polarised proton- 
proton scattering. Such experiments would be possible at RHIC in 
the near future, details can be found in ref. ||9| . 

Actually there are many interesting quantities which can be deter- 
mined in such experiments, also for transverse polarisation, but their 
theoretical understanding is at least partially incomplete and the re- 
liable simulation of background reactions will require fundamentally 
modified codes. We therefore concentrated on the measurement of lS.g 
via direct photon assymmetries, which is the most straight forward 
experiment proposed so far [10|. 

Still it is generally accepted (and will be demonstrated in this con- 
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tribution) that planning and analysing such experiments requires a 
full-fledged Monte-Carlo code. We developped one such code to de- 
scribe the collision of longitudinally polarised nucleons. It is called 
SPHINX and is basically a polarised version of PYTHIA |11|. A 
description of this code will be published elsewhere |12|, and we shall 
try to make it generally accessible. In our paper we present results 
obtained with this code for experiments planned at RHIC using the 
STAR and/or PHENIX detector. 

We present our results in three sections. In section 2 we discuss 
prompt-7 production per se. In section 3 we analyse the additional 
information which is obtained by simultaneously detecting the away- 
side-jet, and in section 4 we present the photon asymmetry generated 
by an asymmetry in the production of vr'^'s subsequently decaying into 
photons. 



2 Prompt-7-Production 

The goal of the measurement of prompt-7-production at RHIC is the 
determination of the gluon polarisation Ag. For this it has to be clar- 
ified that the signal is indeed proportional to /S.g and can be clearly 
separated from the background. For this purpose we investigated the 
two leading processes (i.e. first order in a^) for prompt-7-production, 
namely the Compton process (see figure 1) and the annihilation pro- 
cess (figure 2) and determined their contribution at RHIC energies to 
the cross section for different parton parametrisations. Actually the 
hard matrix elements for prompt-7 production have been calculated 
to NLO but PYTHIA is set up in such a way that the higher 
orders are effectively taken care of by the initial and final state show- 
ering. As the NLO amplitudes do not show any features qualitatively 
different from those of the LO ones this procedure should be fine. 
Note however, that for heavy-quark production the spin-effects could 



change substantially between LO and NLO [14|, which would require 
a more carefull treatment. 

Furthermore we had a closer look at the main contributions to the 
background and analysed the procedures proposed in Q to discrimi- 
nate it. 
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2.1 The Compton process vs. the annihilation 
process 

The hadronic cross section for prompt-7-production in proton-proton 
collisions is given by a convolution of the parton distributions and the 
partonic cross section. In the spin averaged case (o' = 5(cr(|j) + cr(||)), 
where 11 (II) denotes antiparallel (parallel) spins of the two protons) 
it reads 

E,^^f^{s^x,,p^) (1) 

= X] / dXb PaiXa,Q^)Pb{xb,Q'^)Ej ^'^^l~*'^'^ {s,Xf,P±), 

ab 

while for the spin difference (Aa = i (c7(|i) - fT(||))) it is given by 

E,^^^^^is,XF,p^) (2) 

= X/ dXadx,APa{Xa,Q^)AP,{xb,Q^)E^^^^^^{s,XF,P±). 
ab ^ 

Here the sum is over all partonic subprocesses which contributes to 
the reaction pp 7X. Pa and Pf, denotes the unpolarised parton 
distributions of quarks and gluons, AP^ and AP^ the polarised ones. 
The latter are the difference between partons of the same helicity as 
the hadron and those of opposite helicity, the former are the sum of 
the two hclicitics. The partonic cross sections in the helicity averaged 

case is defined as a = ^ {^++ + ^^ ) where the indices +, — signifies 

the helicities of the incoming partons. For Ao" the plus sign between 
the two terms in the formula above has to be replaced by a minus sign. 
Xa and Xf, arc the Bj0rkcn x- variables of the incoming partons, p^ is 
the transverse momentum of the outgoing photon in the pp-CMS. The 
longitudinal momentum fraction of the 7 is defined as = 2p^/i/i. 
The polarised partonic cross section is given by 



dAdab^lX 1 

^7 7^ = otag- 

d-^p-y s 



AM, 



ab^-yX 



(5 (s + (3) 



To obtain the unpolarised cross section one has simply to replace 
the polarised matrix element AMab-^jX by the unpolarised Mab-^'yX ■ 
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Finally the partonic Mandelstam variables are related to the usual 
hadronic ones by: s = XaXbS, t = Xat, and u = XbU. 

In leading order perturbation theory in aog only the Compton 
process qg — > 'yq and the annihilation process qq — > 75 contribute 
to the prompt-7-production. Their polarised and unpolarised matrix 
elements are summarised in table 1. The contribution of the two pro- 
cesses to the hadronic cross section corresponding to (|l]) resp. @ de- 
pends strongly on the parton distributions. In unpolarised pp-collisions 
~ in contrary to pp-collisions - the Compton process clearly dominates 
over the annihilation process, because the gluon density in protons is 
much higher than the antiquark density {g ^ q). In the polarised 
case, however, the relative importance of the two processes depends 
crucially on the relative size of the polarised gluon distribution Ag and 
the polarised sea distribution Aq. Both of them are presently com- 
pletely unknown. By the time the RHIC-Spin-Collaboration (RSC) 
could possibly start to take data, A^(x) should, however, be known 
quite accurately from semi-inclusive lepton-nucleon scattering experi- 
ments by the HERMES collaboration (HERA) |l5|. This should allow 
to avoid ambiguities in the interpretation of potential RSC data. 
For our simulation we used two parametrisations for parton densi- 
ties with large gluon polarisation by Altarelli&Stirling [16] and by 



Ross&Roberts (set D) |17] and one parametrisation with a small gluon 
polarisation by Ross&Roberts (set A). 

For large gluon polarisation the Compton process is the by far domi- 
nant one and one can safely neglect the contribution of the annihilation 
process in (§). In this case the prompt-7-production becomes propor- 
tional to Ag and is such a clean probe for the gluon polarisation: 

dAopp^y^ 
— ¥p^f — 

^ T.J dxa dxb(^Aq{xa,Q^)Ag{xt,Q^)E^^^^§^is,XF,p^) + (x 

However, in a scenario with a large sea contribution to the spin of the 
proton and a gluon polarisation only due to Altarelli-Parisi evolution, 
as described by the parametrisation Ross&Roberts set A [17|, the 
annihilation process becomes the major contribution. 

We investigated these different scenarios with the Monte-Carlo pro- 
gram Sphinx, which can be used to simulate longitudinal polarised 



6 



pp-scattering. We generated 10'' events for both spin combinations 
of the protons at the RHIC energy y/s = 200 GeV. For the unpo- 
larised parton distributions we haven chosen the parametrisation of 
Gliick, Reya, and Vogt |21|, while for the polarised distributions we 
used the parametrisations mentioned above, namely Altarelli&Stirling 
and Ross&Roberts set A and set D. In Sphinx matrix elements are 
implemented in leading order only. However, due to the initial and 
final state shower algorithm some features of higher order effects are 
incorporated as well |11]. Also the polarisation effects are traced in 
the initial state shower. For the simulations the polarised initial state 
shower and the final state shower were switched on. To avoid in- 
frared divergences the hard interaction cross section must be supple- 
mented by a lower cut off for the transverse momentum p± . We chose 
> 4 GeV. 

The results of these simulations are shown in figures 3 to 21. In 
figure ^ the Lorentz-invariant cross section for prompt-7-production 
as a function of p_L at rrj? is displayed for the Compton process 
(upper plot) and the annihilation process (lower plot). is here 
the longitudinal momentum fraction of the photon defined hy xp = 
2pl/y/s. In both cases the spin averaged cross section (squares) and 
the cross section for the spin difference (triangles) are shown. For the 
latter the parametrisation of Altarelli&Stirling has been used. For the 
annihilation process — Ac is plotted, because Act is negative, meaning 
that the cross-section for antiparallel spins is smaller than for parallel 
spin. This can be seen from table 1, keeping in mind that t = xix^t = 
— s(l — cos0)/2 is negative. This negative polarised partonic cross 
section is than multiplied by the positve polarised quark and antiquark 
distributions in (^). The error bars reflect the MC error. A typical 
RHIC run has 320 pb~^, such that the 10^ events we generated for each 
spin combination (with p_L > 4 GeV) corresponds to an integrated 
cross section of 3 • 10^^ mb for the spin averaged case respectively to a 
differential cross-section of roughly 3 • 10~^ mb/(47r • 4 GeV) = 6 • 10^^ 
mb/GeV in the 4 GeV bin. This implies that the MG error is roughly 
comparable to the expected experimental error for the prompt-7's 
and substantially larger than the anticipated experimental errors for 
gammas from vr^ decays. 

At this point we want to state as clearly as possible that we do 
not attribute special significance to any of the used parametrisations. 
In fact virtually nothing is known about Ag{x, Q^), except the trivial 
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fact that its absolute magnitude is limited by g{x,Q^). It could e.g. 
very well be that Ag{x, Q^) changes sign for some x value and this was 
actually advocated to get a good fit to the data in specific models. It is 
just because nothing is known about Ag{x, Q^) that it is so important 
to measure (and RHIC seems to be the only facility able to do so). 
In this situation all our simulations are meant as illustrations. We 
used two similar models (Altarelli&Stirling and Ross&Roberts D) with 
optimistically large gluon-distributions, because it is easy to estimate 
the results for smaller Ag by just scaling them down and for large 
Ag the MC simulation requires less statistics. At the same time the 
difference between the predictions of these two models gives a feeling 
for the sensitivity of experimental signals to details of the distribution 
functions. 

Comparing the two processes one can clearly see that in the spin 
averaged case the Compton process leads to a cross section by an order 
of magnitude larger than the annihilation process. This domination 
is even more pronounced in the polarised case, where one and a half 
orders of magnitude are between both processes. However, the latter 
fact is, as said above, due to the choice of a parametrisation with a 
large Ag and a small Aq. 

In figure |^ we examine the polarised case in more detail by using 
different sets of polarised parton distributions. In the upper part the 
cross section of the spin difference for the Compton process for three 
different sets of parametrisations is shown. The parametrisations of 
Altarelli&Stirling (triangles) and Ross&Roberts set D (squares) have 
a similar, large Ag, but the x-dependence Ag(x) is different, which 
leads to much bigger cross section (logarithmic scale!) in the latter 
case in comparison with the former. This means that it is possible 
to distinguish between both parametrisations in the experiment. The 
third parametrisation Ross&Roberts set A (circles) has a gluon polari- 
sation only due to Altarelli-Parisi evolution and leads to no measurable 
signal. 

In the lower part the annihilation process is shown. Here the results 
for Altarelli&Stirling (triangles), as a representative of the large po- 
larised gluons and small polarised sea scenario, and Ross&Roberts set 
A (circles), as a representative of the opposite situation of nearly un- 
polarised gluons and large, negative polarised sea, are given. As seen 
above, Altarelli&Stirling leads to a small, negative signal, whereas now 
Ross&Roberts set A produce a relatively large, positive signal, which 
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exceeds the corresponding result for the Compton process by far. On 
the other hand it is still by roughly a factor of five smaller than the 
result of Compton process in case of a large gluon polarisation. 

In figure ^ the resulting asymmetries Acr/cr are shown. For the 
Compton process in the upper part the asymmetry grows with for 
Altarelli&Stirling and Ross&Roberts set D and reaches values around 
20% resp. 30%. For Ross&Roberts set A it is consistent with zero. For 
the annihilation process in the lower part Altarelli&Stirling leads to 
small, negative, Ross&Roberts set A to a large, positive asymmetry. 

The main conclusions one can draw at this point are that in the 
case of a large gluon polarisation the annihilation process is a small 
correction, which is calculable using the data from electron-scattering, 
such that prompt-7-production is a clean probe of lS.g. If the absolute 
value of A(7(x) would be a factor ten smaller than assumed by e.g. 
Altarelli&Stirling, which would be the case if the total spin carried by 
gluons were less than half a unit of then annihilation and Compton 
graph would contribute at the same level. It would still be possible to 
extract /^g{x, Q"^), but only by a combined fit to the data from various 
spin-experiments. If A^f would be substantially smaller than 0.5 h the 
determination of /^g{x, Q^) were probably very difficult (unless rather 
large positive and negative parts of lS.g[x) canceled to give a small 
Ag). However, deriving such a low bound for ^g{x) would be very 
interesting as one would not expect it to be that small. The crucial 
remaining questions are whether the anticipated statistics is sufficent 
to actually determine A(^(x) and whether there are background effects 
which could blur the simple picture. These questions will be addressed 
next. 

To examine the precission of prompt-7 measurements at RHIC, 
we transformed these cross sections in counting rates at RHIC and 
determined the statistical errors. To obtain the total rates one has to 
multiply the cross section with the integrated luminosity for which we 
have taken the design value from [^]: / £dr = 3.2 x 10^^ cm~^, which 
corresponds to a luminosity of £ = 8 x 10^^ cm~^ s^^ and an effective 
run time of r = 4 x 10® s, which means 100 days with 50% efficiency. 
At RHIC the beams are only partially polarised: Pbeam = 0.7, whereas 
Sphinx simulates fully polarised events. Therefore one has to combine 
the MC-rates to the experimental rates as follows {P = -Pbeam X-Pbeam): 
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Nil,. = ^A^2/c + ^^<'c- (6) 



For the asymmetry follows: 



]\[U — jsfU 



-iVexp. I -'•exp. 
= i'^MC, (9) 



with the error: 



Mexp ^ ^ ' (10) 

iVJip. + iVjJp. ^Nli. + Nil, 



N^p. + iVjJp. 



exp. 

(11) 



At RHIC several cuts on the events have to be applied. One cut, 
which we have investigated further, is due to the finite coverage in 
the plane of the pseudorapidity r] and the azimuth (p of the detectors 
STAR and PHENIX g. STAR has a full coverage in the azimuth: 
A0 = 27r and a coverage in the pseudorapidity of \rj\ < 1 without end 
caps and of |r/| < 2 in the extended version with end caps. PHENIX 
covers only the half azimuth A(f) = it and has the same 77-coverage 
as STAR. To examine how severe these cuts are we determined the 
rates for prompt-7-production assuming a full azimuthal coverag^ 
and applying the cuts \rj'^\ < 1, \r]'^\ < 2, and no ry'''-cut. The results 
are given in the figures 1,0. The error bars reflect the statistical errors 
at RHIC. 

In figure ^ the spin averaged rates (upper left plot), the spin differ- 
ence rates (upper right plot), and the asymmetry (lower plot) are dis- 
played for the three ?7'^-cuts. The parametrisation of Altarelli&Stirling 
has been used. One can clearly see that the cut \r]'^\ < 2 is not very se- 
vere, especially for high p±, whereas the cut \ri'^\ < 1 reduces the rates 
roughly by a factor of two. In addition also the asymmetry decreases 



-'^This means that we are simulating the situation for STAR. To obtain the corresponding 
results for PHENIX one has to divide the rates by 2. 
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in this case slightly. However, even with the stronger cut Ir^'^l < 1 
similar parametrisations as Altarelli&Stirling and Ross&Roberts set 
D can be distinguished within the experimental errors at RHIC, as 
can be seen in figure |^ Here the rates for the spin difference and the 
asymmetries for the different parametrisations in dependence of the 
applied cuts are compared. Hence at this point the extended version 
of detectors with end caps is not absolutely needed, but, as shown 
later, this extension will be crucial when measuring the prompt-7 and 
the away-side jet. 

2.2 Background considerations 

High-p_|_ 7's are not only produced in the direct processes discussed in 
the subsection above, but at a far larger rate due to bremsstrahlung 
and in particular in meson decays. This background has to be sep- 
arated from the direct photons very accurately in order to do not 
contaminate the signal substantially. In this publication we concen- 
trate on the background produced in pion and r] decays vr'' — > 27 resp. 
r] ^ 2j which is the major contribution and analyse the ability to 
remove it with STAR and PHENIX. 

The mesons which give rise to the background are produced in all 
QCD parton processes, among them the QCD-Compton process: 

qg ^ 9Q — ' . . . 7r° ^ . . . 77, 

qg ^ gq — ' ... ^ ... 77 

is the most important, followed by qq — > qq and gg gg. Beside 
the QCD-Compton process we had a closer look at the i^^-scattering 
process, because this process is especially interesting with regard to 
the gluon polarisation. Here Ag enters twice in the cross section, and 
hence this process is extremly sensitive on changes in Ag. For our 
studies we simulated 10^ events for each polarisation combination of 
the protons for both, the QCD-Compton process qg — > gq and the 
(7g'-scattering process gg gg. 

Comparing the yields of the QCD-Compton process (full symbols) 
with the prompt-7 Compton process (open symbols) in figure |8| it is 
obvious that the background is very important in the case of a large 
gluon polarisation (Altarelli&Stirling left, Ross&Roberts right). In 
the spin averaged case (upper plots) it is much larger than the sig- 
nal, but it decreases also faster with pj_. In addition, the background 
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shows also an asymmetry and the resulting rates for the spin-difference 
(lower plots) are higher than the true signal up to a transverse mo- 
mentum of pj_ ~ 15 GeV. Because the QCD-Compton process is just 
one contribution, although the most important, to the background, 
it is obvious that the photons from neutral meson decay have to be 
separated very accurately from the direct ones in order not to con- 
taminate the true signal. Doing so has the additional advantage that 
the asymmetry in the photons from e.g. pion decay can be used as an 
excellent signal for Ag{x) (see section 4). 

To investigate the background further let us consider how it is com- 
posed. Figure |9| and IC show the compositon for the QCD-Compton 
process and the g'^-scattering process. In both cases the the pion decay 
TT^ —>■ 27 give rise to the main contribution of 80%-90% (logarithmic 
scale!) over the whole p_L-region for both, the spinaveraged case (up- 
per plots) and the spin difference (lower plots). The next important 
contribution is the r/-decay — > 27 with 10%-20%. 

Comparing figure 9 and figure 10 one realizes that the rate of 
the QCD-Compton process is several times larger than that of the 
(73-scattering process in the spin averaged case (upper plots). For the 
spin-difference they are comparable at low p± even with the parametri- 
sation of Altarelli&Stirling with its large gluon polarisation. Due to 
the much faster decrease with p± of the latter the QCD-Compton pro- 
cess becomes also in the polarised case the dominant one at high-p^. 
This dominance will be more pronounced for parametrisations with a 
smaller Ag. 

From the discussion above it follows that the mesons from which 
those background photons stem have to be reconstructed, in order that 
they can be separated from the true signal. There are two main pos- 
sibilities that the mesons escape their reconstruction and thus 'fake' 
photons remain. The two possible sources of 'fake' 7's are the follow- 
ing. First, asymmetric decay, i. e. one 7 is inside the detector, the 
other outside, such that they cannot be combined. Second, merged 
7 's, i. e. the two 7's are too narrow and cannot be resolved by the de- 
tector. Investigations have been done only for the latter case, because 
this is more general, whereas the former is dependent on the precise 
detector geometry and requires a specialized detector simulation. 

In the following we determine the rate of 'fake' 7's from pion decay 
in dependence of the spatial detector resolution. The minimal opening 
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angle of a 7-pair in the rest frame of the pion is given by: 

Xmin = 2^. (12) 

The following resolutions of the detector are considered 

X'""' > 0, 005 rad 
X'"" > 0, 01 rad 

X'"' > 0, 02 rad, (13) 

For the PHENIX detector the planned design would result in x^'^^ > 
0, 01 rad. We define the fake-7-rate R as the fraction of the number 
of unresolved pions and the total number of pions: 

jyunres 

R = ^ (14) 

The MC fake-7-rate Rmc has to be transformed to the experimental 
i?exp as follows: 

oTi(TT) , „oTT(Ti) 



'■exp 

with 



l-P 1-AmC 

P = X : • (16) 

^ 1 + P 1 + ^MC ^ ' 

Then the number of real direct photons is given by: 

= Ntot-RN^. (17) 

The total number of fake-7's in dependence of their opening angle and 
the fake-7-rates for the different resolutions considered of the detector 



are displayed in figure |11|. These plots show that the PHENIX resolu- 
tion of x'^^^ > 0, 01 is sufficient to keep the fake-7-rate below 10% up to 
a transverse momentum of p_L ~ 20 GeV. Due to the steep slope of the 
7-rates around x = 0, 01 in the upper plots the resolution x"^*^*^ > 0, 02 
is absolutely unsuited, whereas the better resolution x'^*^'' > 0, 005 is 
excellent. Further studies are needed to determine the optimal com- 
promise between costs and fake-7-rates for the detector. In addition 
there is a slight spin dependence of the fake-7-rates R^^ < R'^^ such 
that unless x'^'''^ = 0, 005 a carefull simulation is needed to determine 
the resulting corrections and thus keep the systematic uncertainties 
small. 
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3 Prompt-7 and Away-Side Jet 



Additional information can be extracted from prompt-7 measurements 
if they are observed in coincidence with jets, which could also help to 
reduce unwanted background. The resulting information could e.g. be 
used to determine the x-dependence of the gluon polarisation Ag{x) 
which is not sufficiently determined by the prompt-7 signal alone. In 
practice, however, one will determine all the polarised distribution 
functions by simultaneous fits to all data. In the context of such fits 
prompt-7's plus jet data could he of great importance for Ag{x, Q^). 
Let us describe first how we handled the jet-reconstruction. Our 
Monte-Carlo code just like PYTHIA uses the jet routines of Jetset. 
Jet reconstruction was thus done by the jETSET-subroutine LUCELL 
[^ ]. This routine is also used to analyse unpolarised proton-proton 
collisions and defines jets in the two-dimensional {rj — (p)- plane, r] 
being the rapidity and (j) the angle around the z-axes. For our calcula- 
tions we used 25 ry-bins and 24 (p-h'ms with various bounds for \rj\. The 
jet defining algorithm works as follows. First all transverse energies 



E± = \Jp\ + mn? in a bin are summed. If this sum exceeds a certain 
value (E'j^''), which we chose to be 1.5 GeV, than it is treated as a jet 
candidate. Starting from the cell with the highest transverse energy 
all cells in a smaller 'distance' than R = •\/(Ar/p~T~(A^p' (which we 
chose as 0.7) are combined to a 'cluster' and if the total transverse 
energy in this cluster exceeds E^^^^^"^ which we chose to be 3.5 GeV, 
the contained particles are accepted as a jet. 

Finally we use an additional constraint on |r/| implied by the present 
design of the STAR detector. This detector is planned to be built with 
so-called end-caps, giving them a wider r/-range. With these end-caps 
it covers the range \rf\ < 1, r/'^''* < 0.3 while without them only the 



range \rp\ < 2, r]'^'^^ < 1.3. is accessible. 

We start by showing that the larger r] range is really needed to 
do allow for a sensible jet analyses. Figure 12 shows a rapidity distri- 
bution of the jets generated by our code and the consequences of the 
rapidity cuts. In each case the area of one of the small rectangels mea- 
sures the logarithm of the number of events. Obviously the rapidity 
cuts without end-caps are too restrictive while adding the end-caps 
allows to cover nearly all of the interesting rapidity-range. 
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We are interested in the Compton-process of figure 1. Thus we 
require that the jet and the photon are detected with a relative angle 
above 90 degrees in the partonic center of momentum system. In the 
following, jets fullfilling this criterium are called away-side jets. In the 
hadronic center of momentum system, which for RHIC coincides with 
the laboratory system this condition looks rather complicated due to 
the Lorentz-boost 

4x1X2 cos((^T - (^J^*) + (X2e^^ - e-'?'')(x2e'?''* - e-^''*) „ 

COS V = — <1 {) 

^ (X2e^^ + e-"^ ) (X2e^''^* + e-^'^' ) 

(18) 

but it turned out that for the jets we generated this criterium is ac- 
tually equivalent to the much simpler one 

cos{(j)^ -(j)^"') <0 . (19) 

The consequences of this criterium are shown in figure 13. In this 
figure we show histograms of the generated jets as a function of the 
photon and jet rapidity. If the 'away-side-jet' criterium is not used 
we get the results on the left side for the spin averaged and spin 
difference rates. Obviously the photons are in general very strongly 
correlated with the jet axes. If the criterium is applied to single out 
the hard Compton processes the distributions on the right side are 
obtained. Obviously these are only a small fraction of the total events 
and the rapidities are substantially different. The weak remaining 
correlation is a consequence of the Lorentz-boost. In the partonic 
center of momentum system the rapidities are anti-correlated. 

The cross sections for photon-plus-jet events are related to the 
distribution functions according to 

dPajpp ^ 7 + Jet + X) v;- 2 / x / ^ da(gff -yq) 

d'^-(PP-l^ + ^f + ^) ^ Ve^A,(x„)A,(x,)^^^i^^^ 

From this we can calculate the differential 7-jet cross sections as a 
function of if and ry^*^*. Figure 14 and 15 show the results for spin aver- 
age and spin difference in the rf=0, 77''^* =0 bin for our two parametri- 
sations with large lS.g. 

Photon-jet experiments give actually more information than con- 
tained in figure 14 and 15. It is possible to reconstruct the x-values 
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from the measured pseudorapidities and transverse momenta. In the 
ideal case, neglecting all initial and final state interaction this connec- 
tion is simply given by 

Xl (P,0,0,P) 
X2 (P,0,0,-P) 

p_L (cosh rp , cos (jp , sin (jp , sinh rf) 
p± (cosh ri^^^ , cos (f)^^* , sin c^^^^ , sinh r/''^*) . (21) 

where pi and p2 are the four-momenta of the incomming partons. 
Energy-momentum conservation implies than (neglecting all masses) 

Xl ~ 
X2 ^ 

Finally to decide which of these x-values belongs to the gluon and 
which to the quark the following procedure was suggested: Define 
Xa = min(xi,X2) and Xf, = max(xi,X2) and require > 0.2. As the 
gluon distribution is already small at such x-values one can expect 
that .T5 is the quark momentum fraction. 

In proceeding like this a number of rather severe asumptions were 
made such that it was rather unclear how good it would work. SPHINX 
gives us the possibility to check it explicitely. Figure 16 to 19 show 
the results. Figure 16 shows just histograms of the generated quark 
and gluon momentum fractions. Figure 17 shows how these are cor- 
related with Xa, Xb, x^^P, Xg^P. Again the area of the rectangles is 
proportional to the logarithm of the rate, such that the correlations 
are much stronger than they look. The figures on the left side show 
how the X values of the distribution functions x'^^'-^ are correlated to 
those generated by the complete SPHINX algorithm, i.e. the differ- 
ence between x^^^p and Xg respectively x^^p and Xq is entirely due to 
the intrinsic transverse momentum and to effects of the initial and 
final state showering. The identification of the jets is taken from the 
Monte-Carlo, such that there are no misidentifications. In an ideal 
experiment ,t^^p and x^^p are the best measurable approximations for 
Xg and Xq. The right side shows how good the simple procedure just 
described is able to reconstruct the x-values, still without the x^-cut. 



P 



Pi 

P2 

P 

Jet 



2p± 
2p± 



(22) 
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Obviously the quark momentum fraction can be reconstructed quite 
reasonably, while the reconstruction of Xg is problematic. This is why 
a X}, cut is needed. Introducing it the correlation becomes much bet- 
ter, as shown in more detail in figure 18 and 19 for g{x) and Ag{x). 
In these figures the upper graphs shows the histogram of the x values 
actually chosen by the Monte Carlo code. The lower ones show the 
distribution of ar-values reconstructed with the described procedure. 
The agreement is very good for x > 0.04 and gets rapidly bad if one 
goes to smaller a;- values. The problems are far less pronounced for the 
polarised case because the chosen function for Ag{x) is comparably 
small at small x. For a situation where Ag{x) would be concentrated 
at extremely small x it could not be deduced from the 7-jet signal. 
The latter is however true for all polarised experiments. Such a very 
soft polarised gluon distribution would also not be detected in deep in- 
elastic lepton-nucleon scattering such that it could not help to explain 
the observed data. Its only effect could be to screw up the extrap- 
olation to small- a; needed to derive experimental values for the sum 
rules. If the observed data are interpreted as giving evidence for an 
anomalous gluon contribution (however one is trying to define it) than 
Ag{x) cannot be too soft and thus should be seen by RHIC. 

With all these caveats one should note, however, that OTir simula- 
tion shows, that an actual prompt-7-plus-jet experiment would be able 
to distinguish even between the two rather similar gluon distributions 
we used (see figure 19, the bottom plots). 

We conclude this paragraph by stating that prompt-7-jet coinci- 
dences will give interesting data but that their analyses will be highly 
non-trivial requirering extensive numerical simulations. A cut like 
iCfc > 0, 2 is necessary. While a more restrictive cut leads to a better 
reconstruction it also worsens statistic. We hope that SPHINX will 
help to find the best compromise. 

4 7r°-Production 

We discussed already in section 2 that the tt'^ decaying into two pho- 
tons have to be reconstructed in order to extract the direct-7 signal. 
This should be possible with high efficency. However, this opens also 
the possibility to use the 7r° asymmetry as an independent measure- 
ment of the gluon polarisation as well. Their advantage in comparison 
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with the prompt-7's is the larger cross section and hence the smaher 
statistical error. On the other hand the pions are a less direct probe 
for A(7 because several processes contribute to their rate, such that 
the total observed asymmetry depends in a rather involved manner 
on the gluon polarisation. For example the most important contribu- 
tion, the QCD-Compton process is proportional to A^f, but the next 
important one qq qq is independent of Ag, whereas in gg gg it 
enters twice. Again a detailed Monte-Carlo simulation should allow 
to relate the data to the distribution functions. 

In figure the contributions of the QCD-Compton process and 
the g'(7-scattering process to the vr'^-production are displayed. In the 
spin averaged case the QCD-Compton process is the far more impor- 
tant. Due to the large gluon polarisation in Altarelli&Stirling, which 
enters in the gfg-scattering process twice, this process gives the major 
contribution at low p±. However, because it decreases much faster 
with p±, for p_L > 8 GeV the QCD-Compton process becomes the 
dominant one in the polarised case too, even for large gluon polarisa- 
tion. Nevertheless, the i^^-scattering process leads to an asymmetry 
much larger than the QCD-Compton process. Here, the MC-statistic 
gives reliable predictions only up to p± < 15 GeV. The different p±- 
dependence of the individual processes could be used to disentangle 
them. 

Figure 21 shows the yield of pions due to the QCD-Compton pro- 
cess in the spin averaged case (upper plot), for the spin-difference 
(middle) and the resulting asymmetry (lower plot) as a function of p± 
for parametrisations with a large gluon polarisation. While figure 20 
showed results for Xi;' = figure 21 shows the rates integrated over 
xp- In comparison with the prompt-7-datas the two parametrisation 
of Altarelli&Stirling (squares) and Ross&Roberts set D (triangles) are 
much harder to distinguish on the base of these rates. Although this is 
unfortunate in this case it shows that one is really analysing an observ- 
able which tests different properties than the prompt-7- measurements. 

5 Summary 

We have analysed three observables, namely direct prompt-7's, prompt- 
7's in coincidence with an away-side jet, and 7's from vr'^-decay We 
found that each of them should allow to obtain significant results in 
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RHIC spin-physics experiments. These results should allow to settle 
the problem of Ag{x). However, each of these signals has its problems: 
For direct photons the background rate in rather high, for prompt-7- 
jet coincidences the jet reconstruction is non-trivial and has, just as 
for the TT^-decay, to rely very heavily on Monte-Carlo simulations. 
Consequently the existence of several independent codes will be cru- 
cial to relate the observed data in a reliable way to the basic physical 
quantities of interest. We hope that our code will contribute to this 
endeavour. 
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Table 1: The squarred matrix elements for the Compton- 
and annihilation process 
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Figure 1: The Compton graph 
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The annihilation graph 
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Figure 3: Spin-average and spin-difference cross sections for prompt 7 
production. (Note that in the lower plot the negative difference was 
plotted.) 

top: The Compton process bottom: The annihilation process 
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Figure 4: Spin-difTerence cross sections for prompt-7 production for var- 
ious parametrisations. (Note that in the lower plot the negative differ- 
ence was plotted for the Altarelli&Stirling parametrisation.) 
top: The Compton process bottom: The annihilation process 
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Figure 5: The asymmetry for prompt-7 production for various 
parametrisations 

top: The Compton process bottom: The annihilation process 
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Figure 6: Influence of the r/-cut on prompt-7 production 

top left: spin average top right: spin difference bottom: asymmetry 
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Figure 7: Consequences of the ry-cut for two different parametrisations 
top left: spin difference, \r]\ < 1 top right: spin difference, \r]\ < 2 
middle: asymmetry, \ri\ < 1 bottom: asymmetry, \r]\ < 2 
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Figure 8: Comparison of signal and background 
top: spin average left: Altarelli&Stirling 

bottom: spin difference right: Ross&Roberts set D 
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Figure 9: Composition of the background for qg — > gq 
top: spin average bottom: spin difference 
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Figure 10: Composition of the background for gg — ^ gg 
top: spin average bottom: spin difference 
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Figure 11: Fake-7-rates 

1. line: distribution of opening angles for the photon pairs 

2. -4. line: fake-j-rates for different resolutions 
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Figure 12: 7+jet-production for different r\ cuts. The areas of the rect- 
angles are proportional to the logarithm of the counts. 
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Figure 13: Histogramms for the rapidity distribution of jet events. The 
photon rapidity is plotted to the backward-right, the jet-rapidity ry^''^* 
is plotted to the backward- left, both in the range from -2 to 2. Spin 
average and spin difference are compared for all jets and away-side jets, 
top: spin average bottom: spin difference 
left: all jets right: away-side jets only 
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Figure 14: Cross section for prompt-7 production with an away-side jet 
for the Altarelli&Stirling parametrisation 
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Figure 15: Cross section for prompt-7 production with an away-side jet 
for the Ross&Roberts parametrisation 



(xq-Xg) — distribution ; l?7^l<2, l?7'^^'l<1,3 
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Figure 16: Simulated (xg-Xg) distribution for the STAR detector with 
end caps. Xq is plotted to the right, Xg to the left backward, 
top: spin average bottom: spin difference 
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Figure 17: x reconstruction for quarks and gluons, idealy x^'^ should be 
strictly correlated with and Xf, and x^^ with x^^p and Xa. The areas 
of the rectangles are proportional to the logarithm of the counts. 
1.+2. line: quark 3. +4. line: gluons left: xmc vs. Xexp 

1.+3. line: spin average 2. +4. line: spin difference right: x^/b vs. Xg/g 
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Figure 18: Reconstructed spin averaged x^-distributions. MC- 
reconstraction implies that cuts on Xq were applied while for 'exper- 
imental reconstruction' these cuts were imposed on Xb. 
top: MC reconstruction left: AltareUi&: Stirling 

bottom: Exp. reconstruction rechts: Ross&Roberts set D 
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^ X 10' 
7000 



6000 
5000 
4000 
3000 
2000 
1000 




MC — reconstruction 

A difference 

Altarelli&Stirling 



A 



I A 



A 



— X 10- 
7000 



O 



6000 
5000 
4000 
3000 
2000 
1000 










MC- reconstruction 




A 


- 7A 




^ A difference 




A 


I A 


Ross&Roberts 




A Set D 




A 


I A 


A 




A 




A 


I A 








:A 




zr 1 


1 1 1 1 1 1 1 1 1 1 1 



0.05 0.1 0.15 0.2 

MC 



0.05 0.1 0.15 0.2 

MC 
Xn 



^ X 10- 

7000 



6000 
5000 
4000 
3000 
2000 
1000 




exp. reconstruction 

A difference 

A 

▲ 

■^A^ Altorelli&Stirling 
▲ 

A ▲ 

V 

^1 I I I I I I I I I I I I I III 
0.05 0.1 0,15 0,2 

xr^ 



^ X 10- 

7000 



6000 
5000 
4000 
3000 
2000 
1000 




exp. reconstruction 
▲ 



A difference 



Ross&Roberts 
Set D 



0,05 0,1 0,15 0.2 



Figure 19: Same as figure 18 for the spin difference 



7Y°-production ° Vs = 200 GeV, Xp 
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Figure 20: vr'^ production to the hard processes qg —>■ gq and gg gg 
top: spin average middle: spin difference bottom: asymmetry 
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Figure 21: vr'' production: comparison of the results for the 
AltareUi&StirUng and the Ross&Roberts D parametrisation 
top: spin average middle: spin difference bottom: asymmetry 



